Decoherence-free quantum dynamics in circuit QED system 
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We study decoherence in a circuit QED system consisting of a charge qubit and two supercon- 
ducting transmission line resonators (TLRs). We show that in the dispersive regime of the circuit 
QED system one TLR can be used as an auxiliary subsystem to realize decoherence-free quantum 
dynamics of the bipartite target system consisting of the charge qubit and the other TLR condi- 
tioned on the auxiliary TLR initially being a proper number state. Our study gives new insight into 
control and manipulation of decoherence in quantum systems. 
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I. INTRODUCTION 

Decoherence remains a major obstacle to experimen- 
tal realizations of quantum computation and communi- 
cation [H, • As is well known, no system can be com- 
pletely isolated from its environment. Interactions be- 
tween the system and the environment create decoher- 
ence. There are some interesting methods to bypass deco- 
herence in quantum information processing. One of them 
is to encode quantum information into decoherence-free 
subspaces and subsystems (DFSs) 0, 0, S, IE S B 0| • Un- 
der certain conditions, a subspace of a physical system is 
decoupled from its environment such that the dynamics 
within this subspace is purely unitary. The DFS is a set of 
all states which is not Experimental realizations of DFSs 
have been achieved in p hoton systems [HI US G2, El , nu ~ 
clea r sp in systems 0, HH HH , and trapped-ion systems 

nam. 

Advances in circuit QED [19-42] opened new prospects 
in nonclassical state generation and quantum information 
processing in the microwave regime. In the circuit QED, 
superconducting circuits are made to act like artificial 
atoms and a one-dimensional superconducting transmis- 
sion line resonator (TLR) forms a microwave cavity. Un- 
like natural atoms, the properties of artificial atoms made 
from circuits can be designed to taste, and even manip- 
ulated in-situ. Because the qubit contains many atoms, 
the effective dipole moment can be much larger than an 
ordinary alkali atom and a Rydberg atom. This allows 
circuits to couple much more strongly to the cavity. This 
large coherent coupling allows circuits to achieve strong 
coupling even in the presence of the larger decoherence 
present in the solid state environment, then one can ob- 
serve the quantum interactions of matter with single pho- 
tons. Hence, circuit QED can explore new regimes of 
cavity QED. Recently, the circuit QED systems have suc- 
cessfully demonstrated strong coupling between a single 
microwave photon and a qubit [2fJ, the implementation 



of a sing le microwave-photon source in all solid-state sys- 
tem |2l| . as well as single artificial-atom lasi ng J22| and 
interaction between two artificial atoms [HI, [24|. More 
recently, the Lamb shift, two- photon Jaynes-Cummings 
model and controlled symmetry breaking |25[ have also 
been observed experimentally in circuit-QED systems 
[26l . l27l |. These give rise to strong experimental sup- 
ports for on-chip quantum optics and quantum informa- 
tion processing. 

In this paper, we are concerned with decoherence in 
a circuit QED system which includes one SQUID-type 
charge qubit acting as an artificial atom and two super- 
conducting TLRs. We show that in the dispersive regime 
of the circuit QED system one TLR can be used as an 
auxiliary subsystem to realize decoherence-free quantum 
dynamics of the bipartite target system consisting of the 
charge qubit and the other TLR. The paper is organized 
as follows. In Sec. II, we propose the physical model 
under our consideration and present its analytical solu- 
tion in the dispersive regime. In Sec. II, we investigate 
DFS of the circuit-QED system. We show how to realize 
decoherence-free quantum dynamics of the bipartite tar- 
get system. We shall conclude the paper with discussions 
and remarks in the last section. 
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FIG. 1: Schematic setup for the proposed circuit-QED sys- 
tem. The SQUID-based charge qubit is coupled with two 
TLRs, TLRA and TLRB, of lengths L a and Lb, respectively. 
The SQUID is placed at the position of the antinode of the 
quantized voltage of TLRA (i.e., L a /2) and the antinode of 
the quantized current of TLRB (i.e., Lb/4), respectively. 
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II. PHYSICAL MODEL 

Let us illustrate our idea firstly. As shown in Fig.l, 
we consider a circuit-QED system in which an SQUID- 
type charge qubit is coupled to two transmission line res- 
onators, TLRA and TLRB, of lengths L a and Lb, respec- 
tively. The qubit is placed at the position of the antin- 
ode of the quantized voltage of TLRA (i.e., x a — L a /2) 
and the antinode of the quantized current of TLRB (i.e., 
xt, = Lt,/ 4), respectively. It can be controlled by the gate 
voltage, which contains the dc part V g dc and the quantum 
part V a generated by the TLRA, and the biasing flux 
which contains the classical part 4> e and the quantized 
part <!>{, generated by the TLRB threading the SQUID. 

In terms of the annihilation operator a (b) and creation 
operator a t (6 t ) of TLRAJTLRB), the Hamiltonian for 
this system reads as H(| 

H = huJa^a + hu^b + 2Ec(2n g - l)a z - Eja x , (1) 

where u> a and u>b are the microwave frequencies of TLRA 
and TLRB, respectively. The last two terms repre- 
sent the Hamiltonian of the charge qubit. Here a z — 
|1)(1| - |0)(0| and a x = |1)(0| + |0)(1| with 0(1) the 
number of Cooper pairs on the superconducting island. 
Eq = e 2 /2Cs is the charging energy with Cs being the 
total box capacitance. n g is the gate charge number and 
Ej is the Josephson coupling energy given by 



C a V a 



2e 



Ej 



EJ cos I 7T 



(2) 



where C g and C a represent the gate capacitance and 
the coupling capacitance between TLRA and the charge 
qubit. V g dc and V a are the dc gate voltage and the quan- 
tum gate voltage generated by TLRA, respectively. EJ 
is the maximum Josephson coupling energy and $o is 
the flux quanta. The total magnetic flux $ threading 
the dc-SQUID is a sum of two parts $ = $b + $ e with 
$ e being the external classical magnetic flux and $b the 
quantized magnetic flux generated by the quantized cur- 
rent in TLRB. 

The quantum gate voltage and the quantized magnetic 
flux associated with TLRA and TLRB can be expressed 
in terms of the annihilation and creation operators of the 
microwave fields in TLRA and TLRB as 



V a = -JpL(a + a1), 



(3) 



where c and I are the capacitance and inductance per unit 
length for TLRA and TLRB, respectively, S is the area 
of the loop of the SQUID, and d the distance between 
TLRB and the SQUID and /io the vacuum permeability. 
Substituting Eqs. © and © into Eq. Q} we get 

H = Huj a a ji a + huJb^b + 2Ec(2n d g c - \)a z (4) 
-hg a (a + o))a z - EJ cos[0 e + (j) b {b + tf)]a x , 



where n dc — C g V dc /(2e) and we have introduced the 



- g - ^a v g 
coupling constant g a , two parameters <pb and 

are defined by 



They 



9a = 2E c C a \/hu!a/{L a c)/(he), 
4>b = noS^hu>b/(L b l)/(2d<i> ), 

4> e = 7T$ e /$ - 



(5) 



For the simplicity, we choose the classical biasing mag- 
netic flux <& e = and work at the charge degeneracy 
point n dc — 1/2. After making a rotation of —ir/2 around 
the y axis, we get the following effective Hamiltonian 



H = hcj a a'a + TkJbb'b — Tig a (a + a')a x 
+EJcos[Mb + ^)]<J z , 



(6) 



which indicates that under the condition <f>b <C 1, we can 
obtain the following approximation Hamiltonian 



H 



= fko a a'a 
E 



i(6 2 



huj b tfb + EJ[l - (/>l(l + 2tfb)/2]a z 
tf 2 )<T z -hg a {a + a))<T x . (7) 



In order to further simplify the above Hamiltonian, we 
change Hamiltonian ([7]) to the interaction picture with 
respect to the free Hamiltonian of TLRB. After discard- 
ing rapidly oscillating terms, the resulting Hamiltonian 
can be expressed as 

Hj = Tiuj a a)a + ^ L a z - hg a (a + a r )a X: (8) 

where the effective energy separation of the qubit is de- 
pendent of the number operator of TLRB = b^b, and 
it is given by the following expression 



cj q = 2EJ[l-<Pl{l + 2n b )/2]/h, 



(9) 



We consider the case of uj q + uj a ^> uj q — uj a , g a . Then 
under the rotating-wave approximation Hamiltonian ([7]) 
becomes 

Hj = fiuj a a) a H — ~2~ az ~ ^ L 9a( a(7 + + (T - a ^)j (10) 

which is a generalized Jaynes-Cummings model which 
describes the interaction between TLRA and the charge 
qubit with the effective energy separation depending on 
the number operator of TLRB. The quantum electric cir- 
cuit of Fig.l is therefore mapped to the problem of a 
two-level artificial atom inside a cavity. 

We study the dispersive regime of the circuit QED, 
where the cavity and the qubit are out of resonance, 
and the qubit-cavity detuning is larger than the coupling 
strength, i.e., A = u> q — oj a 3> g a . In the dispersive 
regime, uj q /uj a < 1 from Hamiltonian © we can obtain 
the following effective Hamiltonian 



H 



fiu)„a'a ■ 



9l 



2 ' LO a 

x[a z a f a+(<7 z + l)/2], 



1 



(11) 



3 



which can be expressed as the following form in the in- 
teraction picture with respect to the first two terms of 
the Hamiltonian 



Ha 



H 



(12) 



where Hq and H± are defined by 
Hq = huj' a a ! a — hxa' i ab^b, 

H x = -fiu' a {a*a+\) + hxb ] b + hxcJatfb, (13) 

where we have introduced the effective frequency for the 
TLA and the cross-Kerr coupling strength between TLA 
and TLB defined by 



9l 



x 

2' 



X 



hut 



(14) 



which indicate that by a careful choice of the parame- 
ters, it is possible to obtain considerable cross-Kerr non- 
linearity. According to recent experimental data in Ref. 
[3(J, E c /h = 5GHz, Ef/h = 8GHz, oj a = 2vr x 6GHz, 
L a = 25mm, C a = c = 0.13fF, and (j) b = 0.1, we 
find the resulting cross-Kerr coupling strength to be 
X = 360MHz. In circuit QED, the lifetime of the charge 
qubit and the transmission line cavity [30] are about 2/is 
and 160ns, respectively. In the lifetime of the transmis- 
sion line resonator, t = 160ns, we can reach a cross-phase 
shift <j> = x T — 9-17 x 2ir. This means that in the life- 
time of the involved subsystems, we can obtain a large 
cross-phase shift between two microwave fields in the two 
transmission line resonators. 

Obviously, Hamiltonian ([13")) is diagonal in the 



(a z , a^a, b^b) representation with the following basis 



1 



to, n) 



(15) 



where n and to take non-negative integers, i = 0, 1, and 
\i) are eigenstates of the qubit operator <j z with cr z |0) = 
-|0) and a z \l) = |1). The states defined by Eq. ^ are 
eigenstates of the Hamiltonian (fl~3|) . i.e., Hs\m,n,i) = 



E n mi\Tn,n,i) with the eigenvalues 



EmnO = H^'a ~ X n ) m , 
E mn l = -h((jj' a - X n )( m ' 



where the second term is the Hamiltonian of the reser- 
voir, the third one represents the interaction between the 
system and the reservoir with a coupling constant Ck , and 
the last one is a renormalization term [45j . Obviously, the 
interaction term in Eq. (|17[) commutes with the Hamilto- 
nian of the system, this means that there is no energy 
exchange between the system and its environment, so 
that the decoherence described by Hamiltonian Q17[) is 
the phase decoherence. A previous work [43| has shown 
that a nonlinear extension of the model Hamiltonian (fl~7|) 
can well describe the phase decoherence in trapped-ion 
systems. 

The Hamiltonian (fl~7|) can be exactly solved by making 
use of the following unitary transformation 



U = exp 



HsT^^bi-h) 



(18) 



Corresponding to the Hamiltonian (TIT)) , the total den- 
sity operator of the system plus reservoir can be ex- 
pressed as 



Pr(t) 



-iH s t 



u- 



-itY c 



xp T (0)[/" 1 e' 



U 

k Ue r " s K 



(19) 



In the derivation of the above solution, we have used 
p T {t) = U- l p' T (t)U, where p' T = e"^* p' T (0)e i6 ^ with 
H' T = UHtU' 1 and p' T (0) = Up T ^)U- 1 , where p T (0) 
the initial total density operator. 

We assume that the system and reservoir are ini- 
tially in thermal equilibrium and uncorrelated, so that 
Pt(0) = p(0) <g> pr, where p(0) is the initial density op- 
erator of the system, and pr the density operator of the 
reservoir, which can be written as pr — IlfcP/c(0) with 
/5fc(0) is the density operator of the fc-th harmonic oscil- 
lator in thermal equilibrium. After taking the trace over 
the reservoir, from Eq. ([T5)l we can get the reduced den- 
sity operator of the system, denoted by p(t) = trRpx{t), 
its matrix elements are explicitly written as 



P(rn'n'i f )(mni) (^) 



P(m f n f i'){mni) 

(0)R ( 
x exp{— i[E-, 



;!i'n'i')(mm) (p) 

- E mni ]t} , (20) 



III. DECOHERENCE-FREE SUBSPACE 



We now consider the decoherence dynamics of the cir- 
cuit QED system under our consideration. We use a 
reservoir consisting of an infinite set of harmonic oscil- 
lators to model environment of the circuit QED system 
and assume the total Hamiltonian [43|, [44| to be 

H T = H s +J2^kb k h+HsJ2 c k(Pl + h) 



1). (16) where R ,„ 

given by 

R 



) (t) is a reservoir-dependent quantity 



im' n' i'){mni) 



(t) = \{Tr R \v- l {E n 



,)e 



-itu k b+b k 



xVk(E m , n >i>)p k (0)V l [E mni ) 



xc 



itu) k b~tb k 



Vk(Emni)}, (21) 



with Vk(x) = exp[a;c fc (6^ - bk)/uj k ]- 

After somewhat lengthy but straightforward calcula- 
tion, we find that R( m 'n'i')(mni)(t) can be factorized as 



(17) 



■^{m 1 n' i')(mni) \p) 



exp [-i^( m '„'i')( mn i)(t)] 

x exp [-r( m /„/ i /) (m „ i )(f)] . (22) 
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Here the phase shift and the damping factor are defined 
by 

S(f>(m'n'i')(mm){t) = ( E m'n>i' ~ E mni) Ql (*)> ( 23 ) 
^ (m'n'i')(mni) (U = {E m 'n'i' E mn i) Q2{t) 7 (24) 

where the two reservoir-dependent functions are given by 

Qi(t) = [ dujJ(u)^-^-sin(ujt), (25) 
Jo 

Qa(t) = 2^^JH^Msin 2 ^coth^^6) 

where J(uj) is the spectral density of the reservoir, and 
/3 = 1/kgT with ks and T being the Boltzmann constant 
and temperature, respectively. 

Eqs. (f20|) and ([22)1 indicate that the interaction between 
the system and its environment induces a phase shift 
(t) and a decaying factor T( m 'n'i')( m m){t) 
in the reduced density operator of the system. All neces- 
sary information about the effects of the environment is 
contained in the spectral density of the reservoir [SB], [4(| . 

Eqs. (f2"3"|) and (|2~4")) we can see that the phase shift and 
the decaying factor induced by the environment are de- 
termined by the energy differences of the system under 
our consideration E m > n >i> —E mn i and the spectral density 
of the environment J {(J). This implies that decoherence 
can be controlled and manipulated by properly choosing 
the energy differences of the system. In fact, From Eq. 
(fT6|) we can see that there exist three types of energy 
differences given by 

E m 'n>i ~ E mnO = -hu' a (n + n' + 1) + hxm' 

+hx(n'm' + nm), (27) 
E m 'n'o-E mn0 = fwj' a {n' - n) - %x{n'm' - nm) (28) 
Em'n'i - E mnX = -hu' a (n' - n) + %x{m' - m) 

+hx(n' m — nm) , (29) 

which indicates that when m = m' — u' a /x = m , we 
have E mon 'i E man Q — E nlQn 'Q E man Q — E 7non 'i 
E mo ni — 0. Under these conditions from Eq. (24) 
we can find that the damping factor vanishes, i.e., 
F(m n'i')(m ni)(t) = 0. This means that the subspace 
{|mo, n, 1), |mo, n', 0)} with n and n' being arbitrary non- 
negative integers is a decoherence-free subspace of the 
triple-partite circuit QED system under our considera- 
tion. In other words, for our present triple-partite system 



consisting one qubit and two TLRs if TLRA is initially 
prepared the number state |toq), then quantum dynamics 
of the bipartite subsystem consisting of the charge qubit 
and TLRB will be decoherence-free. This means that an 
arbitrary quantum state of the bipartite system would be 
a decoherence-free state conditioned on the auxiliary sub- 
system TLRA initially being the number state |too). In 
this sense, TLRA acts as an auxiliary subsystem which is 
used to control decoherence of the bipartite system con- 
sisting of the charge qubit and TLRB. Hence, we realize 
decoherence-free quantum dynamics of the charge qubit 
and TLRB. 

IV. CONCLUDING REMARKS 

In conclusion, we have studied decoherence in the cir- 
cuit QED system consisting of a charge qubit and two 
superconducting TLRs. Actully, we have proposed a 
scheme to realize decoherence-free quantum dynamics of 
the bipartite consisting of the charge qubit and one super- 
conducting TLR by using another superconducting TLR 
as auxiliary subsystem. In this scheme one TLR and the 
charge qubit constitute the controlled target system while 
the other TLR is the auxiliary subsystem which acts as 
a tool to control the target system. The whole of them 
forms a triple-partite circuit-QED system. It has been 
found that in the dispersive regime of the circuit QED 
system, decoherence-free quantum dynamics of the bi- 
partite target system can be realized when the auxiliary 
TLR subsystem is initially prepared in proper number 
states. This implies that by controlling and manipulat- 
ing the auxiliary subsystem, one can protect quantum 
system against decoherence. This provides fundamental 
insight into the control of decoherence in circuit QED 
systems. It is believed that our present scheme opens an 
new way to engineer decoherence in quantum systems. 
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